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Abstract 

We reformulate in Lagrangian coordinates the two-phase free boundary problem 
for the equations of Magnetohydrodynamics in a infinite slab, which is incompressible, 
viscous and of zero resistivity, as one for the Navier-Stokes equations with a force term 
induced by the fluid flow map. We study the stabilized effect of the magnetic field for 
the linearized equations around the steady-state solution by assuming that the upper 
fluid is heavier than the lower fluid, i.e., the linear Rayleigh- Taylor instability. We 
identity the critical magnetic number \B\ C by a variational problem. For the cases (i) 
the magnetic number B is vertical in 2D or 3D; (ii) B is horizontal in 2D, we prove 
that the linear system is stable when \B\ > \B\ C and is unstable when \B\ < \B\ C . 
Moreover, for \B\ < \B\ C the vertical B stabilizes the low frequency interval while the 
horizontal B stabilizes the high frequency interval, and the growth rate of growing 
modes is bounded. 

Mathematics Subject Classification (2000). 76E25, 76E17, 76W05, 35Q35. 
Keywords. Rayleigh- Taylor instability, MHD, free boundary problem, critical num- 
ber, variational method. 

1 Formulation 

1.1 Formulation in Eulerian coordinates 

We consider the two-phase free boundary problem for the equations of Magnetohydrodynam- 
ics (MHD) within the infinite slab = x (—1, 1) c M. d , d = 2 or 3 is the dimension, 
and for time t > 0. The fluids are separated by a moving free interface S(t) that extends 
to infinity in every horizontal direction. The interface divides Q into two time- dependent, 
disjoint, open subsets fi±(t) so that Q = Q+(t) Ufi_(t) U£(t) and £(t) = Q+(t)nfi_(t). The 
motions of the fluids are driven by the constant gravitational field along e^— the Xd direction, 
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G = (0, 0, — g) with g > and the Lorentz force induced by the magnetic fields. The two 
fluids are described by their velocity, pressure and magnetic field functions, which are given 
for each t > by, respectively, 

{u±,p±,h±){t,-) : Ot(t) ->• (R + ,R d ,R,R d ). (1.1) 

We shall assume that at a given time t > these functions have well-defined traces onto 
S(t). 

The fluids under consideration are incompressible, viscous and of zero resistivity, hence for 
t > and (x',Xd) <E Q±(t) the fluids satisfy the following magnetohydrodynamic equations: 



1.2) 



d t {p±u±) + div(p ± M ± <g> u±) + divS± = -gp±e d , 
divu± = 0, 

d t h± + div(-u ± <g) h±) — div(h± <g) -u ± ) = 0, 
di\h± = 0, 

where we define the stress tensor consisting of fluid part and magnetic part by 

\h I 2 

S± = -/i±(Vu ± + Vul) + p±I + -h±® h±. (1.3) 

Hereafter the superscribe T means the transposition and / is the d x d identity matrix. The 
positive constants p± and /i± denote the densities and viscosities of the respective fluids. 

Now we prescribe the boundary conditions at the fixed boundaries and the jump con- 
ditions at the free interface. First for the the fluid equations (1.2) 1; due to the viscosity 
we assume that the velocity is continuous across the free interface and enforce the no-slip 
condition at the fixed upper and lower boundaries. Since we do not take into account the 
surface tension, it is standard to assume that the normal stress is continuous across the free 
interface (cf. PQ[H]). Therefore, we impose the boundary conditions at the fixed boundaries 

u + (t,x',-l) =«_(t,x', 1) =0, forallt>0, x'eR^ 1 , (1.4) 

and impose the jump conditions at the free interface 

M |e(*)= 0, (1.5) 
[Su] \ m = 0, (1.6) 

where we have written the normal vector to S(t) as z/, for the trace of a quantity / on 

S(t) and denote the interfacial jump by 

[/] |s(t) : = /+|s(t) - /-|s(t)- (1-7) 

For the magnetic equations (1.2)3, since the fluids are of zero resistivity this is a free transport 
equation along the flow, and hence the Dirichlet boundary condition on the velocity at the 
fixed boundary prevents the necessity of prescribing boundary condition on the magnetic 
fields. On the other hand, due to the divergence-free (1.2) 4 and also physically, we shall 
assume that the normal component of magnetic field is continuous across the free interface 
(cf. Pi) 

[h-u]\ m = 0. (1.8) 



2 



However, we will show explicitly that the divergence-free of h± and the jump condition (1.8) 
can be verified if they hold initially. Therefore, the conditions (1.2) 4 , (1.8) are transformed 
to the compatibility conditions assumed on the initial magnetic field. 

The motion of the free interface is coupled to the evolution equations for the fluids (1.2) 
by requiring that the boundary be adverted with the fluids. More precisely, if V(t,x) G M. d 
denotes the normal velocity of the boundary at x G S(t), then 

V(t, x) = (u(t, x) ■ v(t, x))v(t, x). (1.9) 

Here u(t,x) is the common trace of u±(t,x) onto £(£) and these traces agree because of the 
jump condition (1.5), which also implies that there is no possibility of the fluids slipping 
past each other along £(£). 

To complete the statement of the problem, we must specify initial conditions. We give 
the initial interface E(0) = S , which yields the open sets f2±(0) on which we specify the 
initial data for the velocity and magnetic field 

(u±,h±)(0,-) : Ofc(0) ->• (R d ,R d ). (1.10) 

To simply the equations we introduce the indicator function \ an d denote 

f u = u + xn+ + u^x^i- , h = h + xn+ + h-Xn- , P = P+Xn+ + P-Xn- , 
1 P = P+Xn+ + P-Xsi- , P = P+Xn + + P-Xsi- , 

and also define the modified pressure by 

\h\ 2 

p = p + — + gpx d - 

Hence the equations (1.2) are replaced by 

pd t u + pu ■ Vu — pAu + Vp = h ■ V/i, 
d t h + u ■ Vh - h ■ Vu = 0, 
divu = divh = 0, 

and the jump condition (1.6) becomes, setting [p] = p + — p_, 

[(pi - p{Vu + Vu T )) v] | S (t)= g[p]x d jy + h-u[h] | S (t) • 

1.2 Formulation in Lagrangian coordinates 

Since the movement of the free interface £(£) and the subsequent change of the domains 
fi±(t) in Eulerian coordinates result severe mathematical difficulties, we switch our analysis 
to Lagrangian coordinates so that the interface and the domains stay fixed in time. To this 
end we define the fixed Lagrangian domains Q + = IR d_1 x (0, 1) and f2_ = IR d_1 x (—1, 0). 
We assume that there exist invertible mappings 

ri° ± : Q± ->■ Ofc(0), (1.15) 
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(1.11) 
(1.12) 



(1.13) 



(1.14) 



which are continuous across {x d = 0} so that S = 77° ({x^ = 0}), {xd = 1} = ^({x^ = 1}), 
and {xd = —1} = rf_{{xd = —1})- The first condition means that S is parameterized by 
the either of the mappings i]± restricted to {X2 = 0} (which one is irrelevant since they are 
continuous across the interface), and the latter two conditions mean that i]± map the fixed 
upper and lower boundaries into themselves. Define the flow maps r\± as the solution to 

d t r)±(t, x) = u±(t, r)±(t, x)), ^ ^ 

f]±(0, x) = r]±(x). 

We think of the Eulerian coordinates as (t,y) with y = rj(t,x), whereas we think of La- 
grangian coordinates as the fixed (t,x) G R + x f2. In order to switch back and forth 
from Lagrangian to Eulerian coordinates we assume that that rj±(t, •) are invertible and 
^±(0 = ^±(^^±), and since u± and 77° are all continuous across {x d = 0}, we have 
S(t) = f]±(t,{xd = 0}). In other words, the Eulerian domains of upper and lower fluids 
are the image of Q± under the mappings i]± and the free interface is the image of {xd = 0} 
under the mapping r]±(t, •). 

Setting r] = X+V+ + X-V-i we define the Lagrangian unknowns 

(v,b,q)(t,x) = (u,h,p)(t,r](t,x)), (t,x)eR + xtt. (1.17) 

Defining the matrix A via A T = (Drj) -1 , then in Lagrangian coordinates the evolution 
equations for i], v, b, q are, writing dj = d/d Xj , 

d t r)i = Vi, 

pd t Vi + A jk d k Tij = bjAjkdkb 1 , 

AjAvj = 0, (1.18) 
d t bi = bjA jk d k Vi, 
Ajkdtbj = 0, 

where the stress tensor of fluid part in Lagrangian coordinates, T(v, q), is given by 

Tij = qlij - fi(A jk d k Vi + A ik d k Vj). (1.19) 

Here we have written 1^ for i,j component of the identity matrix / and we have employed 
the Einstein convention of summing over repeated indices. 

To write the jump conditions, for a quantity / = f±, we define the interfacial jump as 

Ul-=f+\{*r4}-f-\lx d =o}- (1-20) 
Then the jump conditions in Lagrangian coordinates are 

H = 0, {bjUjl = 0, [?>,■] = g[p] Vd n l + MbjUj, (1.21) 

where the unit normal to the interface £(£), i.e., n = v{rf) can be represented by 

Ae d 



n 



\Ae rl 



(1.22) 



{xd=0} 

Finally, we require the no-slip boundary condition 

v + (t,x', -1) =v-(t,x',\) = 0, foraUf>0, x' G M d_1 . (1.23) 
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1.3 Reformulation 

One purpose of this paper is to reformulate the free boundary problem (1.18), (1.21), (1.23) 
as the Navier-Stokes equations with a force term induced by the fluid flow map. We want 
to eliminate b by expressing it in terms of n. Indeed, applying An to (1.18)4, we have 

Audth = bjAjkdkViAn = bjA^d^d^An = -bjA^dk^dtAn = -bid t A a . 

This implies that d t (Ajibj) = and hence, 

Ajibj = A%b% (1.24) 
h = d lV l A%b]. (1.25) 

Hereafter, the superscript means the initial value. 

By the expression (1.25), we first check the divergence of b, i.e., (1.18)5. We make use 
of the geometric identities 

J = J° and d k (JA ik ) = 0. (1.26) 
Hence, applying A ik d k to (1.25) we have 

Akd k h = ^AMd^A ^) = ^d k {JA lk dtfA%q) = ^d k {J°A» k b]) = A%d k b». (1.27) 

Next we check the jump [^-n?]. We recall the expression (1.22) of unit normal to the interface 
S(t). It is easy to verify that Aej is continuous across the free interface. Hence we have 

IbjUj} = ld lV >A kl bl ■ A jd }j-^ = |il-7 = ^ n 3^f (1-28) 
Hence, if we assume the compatibility conditions on the initial data 

A%d k tfj = 0, [&X ]=0, (1.29) 
then from (1.27), (1.28), we have 

A jk d k bj = 0, {bjUjl = 0. (1.30) 

Moreover, for simplify of notations, we assume that 

A® m $n = Bi with B is a constant vector. (1-31) 

We remark that the class of the pairs of the data r]°,b that satisfy the constraints (1.29), 
(1.31) is quite large enough. For example, we chose 77° = Id and b° = const, then by (1.24), 
(1.27), (1.28), any pair of data n, b which is transported by the flow will satisfy (1.29), (1.31). 
Now we represent the Lorentz force term by, since (1.24), (1.25), (1.30), (1.31), 

b 3 A 3k d k U = dtrfAlfilA^drfAlb") = A° mk b° m d k (d r ^ A° sr b° s ) = B^d^f. (1.32) 
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Hence the equations (1.18) becomes a Navier-Stokes system with the force term induced by 
the flow map rj: 

dtVi = v t , 

pd tVi + A jk d k Tij - B^d^rt = 0, (1.33) 
AjkdkVj = 0, 

where the magnetic number B can be regarded as a vector parameter. Accordingly, the jump 
conditions (1.21) become 

H = 0, [7>J = g[p} Vd n l + BiBrnldrfjd^rf, (1.34) 

Note that we implicitly admit that B m d m rfr^ is continuous across {xd = 0}, this follows 
from the assumptions (1.29), (1.31). Finally, we require the boundary condition (1.23). 

1.4 Linearization around the steady state 

The system (1.33), (1.34), (1.23) admits the steady solution with v = 0, n = Id, , q = const 
with the interface given by r]({xd = 0}) = {xd = 0} and hence n = e^, A — I. Now we 
linearize the equations (1.33) around such a steady-state solution, the resulting linearized 
equations are 

d t r] = v , 

pd t v + Vq-(xAv-B l B m d? m r ] = 0, (1.35) 
divu = 0. 

The corresponding linearized jump conditions are 

H = 0, l-p(Dv + Dv T ) + qlje d = g[p]v d e d + B d Bi{divl (1-36) 

while the boundary conditions are 

v-(t,x',-l) =v+{t,x',l) = 0. (1.37) 

The main purpose of this paper is to study the stabilized effect of magnetic field on the 
Rayleigh- Taylor problem, hence we assume that the upper fluid is heavier than the lower 
fluid, i.e., 

p + >p_^ [p] > 0. (1.38) 

1.5 Normal mode ansatz 

In the fluid stability analysis it is standard to study the normal mode solutions of the 
linearized system (cf. PQ). To begin, we assume an ansatz 

v(t,x) = w(x)e xt , q(t,x) = q(x)e xt , n(t,x) = f](x)e xt , (1.39) 

for some A > 0. Substituting this ansatz into (1.35), eliminating fj by using the first equation, 
we arrive at the time-invariant system for w = (wi, . . . , Wd) and q: 

Xpw + Vg - pAw - X^BiB^^w = 0, ^ ^ 

divw = 0, 
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with the corresponding jump conditions 

H = 0, l-fi(Dw + Dw T ) + qlje d = \- l g[p]w d e d + \~ l B d Bi {diwj, (1.41) 
and the boundary conditions 

w_(t,x',-l) = w + (t,x',l) = 0. (1.42) 

We make the further structural assumption that the x' dependence of w, q is given as a 
Fourier mode e tx ''^ for £ e R d_1 . Together with the growing mode ansatz, this constitutes a 
"normal mode" ansatz. At the rest of this subsection, we shall write down the analysis for 
the three dimension d — 3, the case for d = 2 can be tracked readily. We define the new 
unknowns </?, 9, ip, ir : (—1, 1) — >■ R by 

Wl (x) = / r -(.r :i ),.' Vi . W2 (x) = -i9(x 3 )e ix '-t, w 3 (x) = ^(x 3 )e tx '<, q(x) = n(x 3 )e ix '<. (1.43) 

To write down the equations for if, 9, ip, n, we first notice that, denoting ' = d/dx 3 , 

(Dw + Dw T )e 3 = - <?'), i(W - 9'), 2ij') T . (1.44) 

The equations will be quite different for the cases that B is vertical and horizontal. 

We first treat the vertical case that B = (0,0, B). In this case, we deduce from the 
equations (1.40) that for each fixed nonzero spatial frequency £ = (£1,^2), and for (p,9,xjj,n 
and A we arrive at the following system of ODEs 



\ 2 fxp - A£ l7 r + //A(|£| V - </) - \B\ V = 0, 
\ 2 P 9 - A£ 2 7r + fi\{\£\ 2 9 - 9") - \B\ 2 9" = 0, 
\ 2 pi) + Avr' + p\{\i\ 2 i) - if/') - \B\ V = 0, 



(1.45) 



along with the jump conditions 



= m = m = 0, 

IpiXfal, - <//) + \B\ V] = - #0 + 1^'] = 0, (1.46) 
[-2^' + A7r-|S|Vl=^, 

and the boundary conditions 

<p(-l) = <p(l) = 0(-l) = 9(1) = V(-l) = V(l) = 0. (1.47) 

Now eliminating 7r from the third equation in (1.45) by using the other equations, we 
arrive at the following fourth-order ODE for ip 

- A 2 P (i£i v - in = A*A(ieiv - 21^1 v + in + i^i 2 (-i£i v + in, a-^) 

along with the jump conditions 

W = M = (1.49) 

[Md^ + O + I^IVHo (1.50) 

[M(^'-3|£|V) + |S|V1 = [AVj + <?[p]|£lV (1.51) 
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and the boundary conditions 



iP(-l) = V(l) = V'(-l) = ^'(1) = 0. (1.52) 

Consequently, for B is vertical to look for growing normal mode solutions of the original 
linearized problem (1.35)-(1.37) reduces to find solutions of the ordinary differential system 
(1.48)-(1.52). Hence, we may say in this paper that if there exists such a pair (A, ip) with 
A > satisfying (1.48)-(1.52), then the linearized problem (1.35)— (1.37) with vertical B is 
unstable; otherwise, the linearized problem is stable. 

Now we treat the horizontal case that B = (5,0,0), without loss of generality. We 
deduce from the equations (1-40) the following system of ODEs for if, 9, ip, 7r and A 



\ 2 p<p - A&tt + M(|C| V - </) + |5| 2 £V = 0, 
\ 2 P 6 - A&tt + /iA(|£| 2 # - 0") + \B\ 2 g0 = 0, 

\ 2 P i) + Att' + Mdei 2 ^ - V") + Wili) = o, 



along with the jump conditions 



(1.53) 



= 0, 

IpXiM - iff)} = [/iA(6^ - e')j = 0, (1.54) 
HW + Att] = 

and the boundary conditions 

<p(-l) = <p(l) = 0(-l) = 6(1) = ^(-1) = V(l) = 0. (1.55) 

Eliminating tc from the third equation in (1.53) we obtain the following ODE for ip 

- A 2 p(iei 2 ^ - r) = Mdei v - 2iei v + in + i^raei 2 ^ - &n (i.56) 

along with the jump conditions 

W = M = 0, (1.57) 
bA(|eiV + ^")]=0, (1-58) 
[/xAty'" - 3|C| V)] = [AW] + </MKlV, (1-59) 

and the boundary conditions 

^(-1) = ^(1) = ^'(-1) = ^'(1) = 0. (1.60) 

2 Main results 

We first remark that if p: — B — 0, i.e. for the two-phase free boundary problem of the 
incompressible Euler equations, the equation for ip is replaced by 

P \^- P r = o, (2.1) 
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along with the modified boundary conditions 

^(-i) = V(i) = o, M = o, W\ = -\- 2 g[ P ]\t\ 2 ^. (2.2) 

We directly obtain from (2.1)-(2.2) that 



1 

'|2\ j„. _ \-2 r ll(N2„/,2/ 



p(|£fV + )&3 = A-^[p]|e|V(0). (2.3) 



i 



This immediately implies that if [p] > 0, then nontrivial solutions ■?/>(£, £3) with A(£) > can 
be found. Moreover, A(£) can be chosen as A(£) > C|£| which implies A(£) growing arbitrarily 
as £ — > 00. This leads to the classical ill-posedness of the Rayleigh- Taylor problem for the 
Euler equations without surface tension, see [21 H] . 

Notice that in the absence of viscosity (of course with the boundary conditions modified 
accordingly) and for any fixed spatial frequency £ 7^ 0, (1.48)-(1.52) (or (1.56)-(1.60)) can 
be viewed as an eigenvalue problem with eigenvalue —A 2 . It is not hard to check that such 
a problem has a natural variational structure that allows for the construction of solutions 
via the direct methods, see below. However, the presence of the viscosity results the appear- 
ance of A both quadratically and linearly which destroys the variational structure of these 
eigenvalue problems. To restore the ability to use variational methods, the authors in [5] 
developed a quite general and robust method. More precisely, they artificially remove the 
linear terms pX by the modified viscosities p = sp where s > is an arbitrary parameter. 
After establishing the existence of the solutions to the modified problem, by proving that 
there is a fixed point A = s, then the corresponding solution is a solution to the original 
problem. We use the same trick in this paper to construct solutions. However, the purpose 
of this paper is to study in detail the stabilized effect of the magnetic field in the Rayleigh- 
Taylor instability, to characterize the critical magnetic number and the critical frequency 
when the magnetic number under the critical value. 

We assume in this paper that p, B are not zero. We define the critical magnetic number 
\B\ C through the following variational problem 

|BE:= sup , (2.4) 

which only depends on g, [p]. For \B\ < \B\ C , we define the critical frequency \C,\^ C for the 
vertical B by the variational problem 

\B\ 2 f\ W'\ 2 dy 
' ^((-1,1)) g[p)^(0) - \B\ 2 f^\iPf dy' 

and define the critical frequency |£|f c for the horizontal B by the variational problem 

(\t\hc) ■■= sup — -3 . (2.6) 

tfeflo 1 ((-1,1)) \B\ 2 J^l^dy 

The superscript emphasizes the dependence of the critical frequencies (2.5)-(2.6) on \B\, 
these two only depend on g, [p] and \B\. We will show in Lemma 3.2 that the extremums 
(2.4)-(2.6) are achieved by the direct method. 
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Now we may first state out main results of the paper for the two-dimensional case, the 
results for three-dimensional case will be discussed later. The first one is concerned with the 
problem (1.48)-(1.52) for the vertical magnetic number. 

Theorem 2.1. Let B = (0,B) be vertical, then we have: 

(i) For any fixed \B\ > \B\ C and any £ G 1R there is no nontrivial solution ip with A > 
to the problem (148)-(1.52); 

(ii) For any fixed \B\ < \B\ C , then for £ G R so that |£| < |£|f c there is no nontrivial 
solution ip with A > to the problem (l.J h 8)-(1.52); 

(in) For any fixed \B\ < \B\ C , then for £ G IR so that |£| > |£|f c there exists ip = ip(^,X2) 
and A(£) > to the problem (1.48)-(1.52). Moreover, ip, A are even in £ and ip is smooth 
when restricted to (— 1, 0) or (0, 1) with 0) 7^ 0. 

Next theorem is concerned with the problem (1.56)-(1.60) for the horizontal magnetic 
number. Note that since we consider the 2D case, £f = |£| 2 for the moment. 

Theorem 2.2. Let B = (B, 0) be horizontal, then we have: 

(i) For any fixed \B\ > \B\ C and any £ G K there is no nontrivial solution ip with A > 
to the problem (1.56)-(1.60). 

(ii) For any fixed \B\ < \B\ C , then for £ G 1R so that \£\ > |^|f c there is no nontrivial 
solution ip with A > to the problem (1.56)-(1.60). 

(in) For any fixed \B\ < \B\ C , then for £ G WL so that |£| < |^|f c there exists ip = ^(^5 ^2) 
and A(£) > to the problem (1.56)-(1.60). Moreover, if), A are even in £ and ip is smooth 
when restricted to (—1, 0) or (0, 1) with ip(£, 0) 7^ 0. 

As we will see, the conclusions of Theorem 2.1 and Theorem 2.2 also hold for the inviscid 
fluids (i.e., fi = 0). They demonstrate that the stabilized effect of the magnetic field is 
much more remarkable than that of the surface tension a which only stabilizes the frequency 

interval (0, y— ) for any a > 0, see [Tj [S]. 




It is easy to imagine the corresponding results for the 3D case. Indeed, for the vertical 
magnetic number B = (0,0,5), the problem (1.48)-(1.52) is the same and hence Theorem 
2.1 holds. However, for the horizontal B = (B, 0, 0), if we consider the frequency £ = (£1, £2) 
with £x = 0. Then (1.56)-(1.60) becomes (1.48)-(1.52) with B = 0, there is no effect of the 
magnetic number and hence by Theorem 2.1 (Hi) the system is unstable. In other words, 
the horizontal magnetic number only stabilizes the frequencies along its own direction. 

It is important to know the behavior of the eigenvalue A(£) obtained in (Hi) of both 
Theorem 2.1 and Theorem 2.2 within their respective definitional intervals. 

Theorem 2.3. Let < \B\ < \B\ C , then we have 

(i) Let A(£) be that of (in) in Theorem 2.1, then A : (|£|^,, 00) — > (0, 00) is continuous 

and 



(ii) Let A(£) be that of (Hi) in Theorem 2.2, then A : (0, |£|f c ) — >• (0, 00) is continuous 

and 





(2.7) 




0. 



(2.8) 
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It is worth to point out that the proof of Theorem 2.3 is independent of the viscosity \i 
and Theorem 2.3 allows us to define the fastest growth rate by, for vertical and horizontal 
magnetic number, respectively, 

A v = sup A(|£|) < oo and A h = sup A(|£|) < oo. (2.9) 

i£i? c <i£i<oo o<iei<i^|B 

This means that the growing mode solutions constructed in Theorem 2.1 and Theorem 2.2 
(except for the horizontal B in 3D case) do not grow arbitrarily, which may imply that the 
magnetic field prevent the ill-posedness for the inviscid Rayleigh- Taylor problem, see [2J H] . 

With if} constructed in (Hi) of Theorem 2.1 (resp. (Hi) of Theorem 2.2), one can use 
the system of ODEs (1.45)-(1.47) (resp. (1.53)-(1.55)) to obtain a growing normal mode 
solution to the linearized problem (1.35)-(1.37) for B is vertical (resp. B is horizontal). 
Although e lx ' ^ L 2 (Q), as in jH |5] , we can resort to a Fourier synthesis of such solution to 
construct solutions to (1.35)-(1.37) which grow in the piecewise Sobolev space H k , Wk > 1. 
The growth rate of such growing solutions is not exactly the fastest growth rate e v (resp. 
e tAh ), but arbitrarily close to this rate. On the other hand, based on the proof of Theorem 
2.1-2.3, we can estimate the growth in time of arbitrary solutions to (1.35)-(1.37) in terms 
of e* A " (resp. e tAh ). The statement and the proof of these results is very similar to that of 
[5j Theorem 2.4-2.5], we omit them and let the interested readers refer to [5] for the details. 

Notice that in Theorem 2.1-2.2, when \B\ > \B\ C no growing mode solutions can be 
constructed by our argument. This should suggest that the linear system (1.35)-(1.37) is 
stable when \B\ > \B\ C . In fact, when \B\ > \B\ C we can prove the following stability 
estimates: 

Theorem 2.4. Let \B\ > \B\ C be fixed, then 

(i) for B = (0, B), we have 

\\d t v(t)\\h + Ht)\\h + \\d*v(t)\\h + f \\dtv(s)\\m ds 
<C(||^(0)||i 2 + |K0)||| 2 + ||^(0)||i 2 ); (2.10) 
if in addition, r](0,Xi, —1) = r](0,xi, 1) = initially, then we have 

\\v(t)\\h+\\v(t)\\l> + \\&m\b + J*\K*)\\*d* 

<C(|K0)||| 2 + 1177(0)111. + 11^(0)1120. (2.H) 

(ii) for B = (B, 0), we have 

\\dtv(t)\\h + ||^i(*)|li> + \\v2(t)\\ 2 m + jf \\d t v(s)\\* H1 ds 

< c(||^(o)||| 2 + \\ehvi(p)\\h + ll^(o)ll^); (2-12) 

if in addition, (dirji + 92r/ 2 )(0) = initially, then we have 

Ht)\\ 2 L2 + Whrnmh + \\m(t)\\h + 1 IK*)II^ da 

<C(\\v(0)\\l 2 + \\d lVl (0)\\ 2 L2 + \\v2(0)\\ 2 m ). (2.13) 
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As before, the assertion (z) also holds for 3D case. The proof of Theorem 2.4 follows by 
combining the standard energy estimates and the definition (2.4) of the critical magnetic 
number \B\ C . It is understood that the initial value d t v(0) are given in terms of the initial 
data through the equations (1.35)2- It is obvious that the estimates (2.10)-(2.12) also hold 
for d\d{, Vijj > 1 since they satisfy the same system (1.35)-(1.37). It in particular implies 
that the exponential growth of solutions is impossible. 

At last, we emphasize that the success of deriving the critical magnetic number (and 
hence critical frequency) essentially depends on the fact that we consider the problem within 
a slab. The critical magnetic number does not hold for the whole space problem and indeed 
in this case for any B there always exists growing normal mode solutions, see pQ. 

The instability of certain steady states of different densities which occurs under an ac- 
celeration of the fluid system in the direction toward the denser fluid is well known as the 
Rayleigh- Taylor instability since the works [HI HQ] . It attracts huge attention of the researches 
both physically and mathematically, we refer to book pQ and the report [7] for the detailed 
references. Our construction of the growing normal mode solutions by variational method is 
inspired by the works [31 SI El E] which deal with the nonconstant density profile. 

The rest of the paper is devoted to prove our theorems. 

3 Proof of Theorems 
3.1 Proof of Theorem 2.1 

We treat the problem (1.48)-(1.52) to prove Theorem 2.1 in this subsection. As quoted in 
Section 2, to restore the ability to use variational methods we remove the linear dependence 
on A by defining the modified viscosities fx = sfi, where s > is an arbitrary parameter. 
Then we obtain a family (s > 0) of modified problems 

- A 2 P (iei 2 ^ - r) = v - 2iei v + n + wc-iei v + <n, (3.1) 

along with the jump conditions 

M = Wl = 0, (3.2) 
MI£IV + V>'0 + |£|V'] = 0, (3.3) 

MiT - 3|e| V) + \b\ Yl = [A 2 W'] + <?[p]|£lV (3-4) 

and the boundary conditions 

^(-1)=^(1) = ^(-1)=^'(1) = 0. (3.5) 

Notice that for any fixed s > and £, (3.1)-(3.5) is a standard eigenvalue problem for 
—A 2 . It allows us to use the variational methods to construct solutions. We define the 
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energies 



em = \ £ s»m 2 m 2 + 1 \^ + r\ 2 ) + i^aei vi 2 + m 2 ) ^ 2 
-\\eg[p]\m\ 2 , (3.6) 



i ^ 



■W) = n / p^iw + ivt)^, (3.7) 
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which are both well-defined on the space Hq((—1, 1)), the subset of if 2 ((— 1, 1)) satisfying 
(3.5). We define the set 

A = {ii>eH*((-l,l)) | JW = 1}, (3.8) 
then we want to find the smallest —A 2 by minimizing 

-A 2 (|el) = «(lel):= jnf/W. (3.9) 

The first thing is to show that a minimizer of (3.9) exists and that the minimizer satisfies 
Euler-Langrange equations equivalent to (3.1)-(3.5). 

Proposition 3.1. (i)For any fixed £ ^ and s > 0, E achieves its infimum on A. 

(ii) Let ip be a minimizer and —A 2 = a := E(ip), then the pair ip, A 2 satisfies (3.1) along 
with the jump and boundary conditions (3.2)-(3.5). Moreover, tp is smooth when restricted 
to (-1,0) or (0,1). 

Proof. To prove (i) we first note that for any ip G A 

m > -\mp]\m\ 2 > -\\e\p\ £ w?^ > -M-^we- (3.10) 

Hence E is bounded below on A. Let ip n G A be a minimizing sequence, then E(ip n ) is 
bounded. This together with (3.6), and (3.10) again imply that ip n is bounded in Hq((—1, 1)). 
So there exists ip G Hq((— 1, 1)) such that ip n ^ ip weakly in ip G Hq((— 1, 1)) and strongly 
in C 1 ((— 1, 1)). ip £ A follows from the strong convergence. Moreover, by the lower semi- 
continuity and the strong convergence, we have 

E(i/>) < liminf E(ip n ) = ME. (3.11) 

n— >oo A 

Now we prove (ii). We notice that since E and J are homogeneous of degree 2, (3.9) is 
equivalent to 

-A*= inf f« (3.12) 

Vei/ 2 ((-i,i)) J(^) 

For r G K and any *0o G Hq{(—\, 1)), define ^(r) = tp + tt/'o, then (3.12) implies 

E(iP(t)) + \ 2 J{^{r)) > 0. (3.13) 
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If we set I(t) = E(ip(r)) + \ 2 J(ip(r)), then we have 7(r) > for all r G R and 1(0) = 0. 
This implies J'(0) = 0. By the expressions (3.6)-(3.7), direct computation leads to 

f s/*(4|e|Wo + (KlV + Vo + <)) + \B?m Wo + «') d*2 

= |e| 2 [p]^(0)^o(0) - X 2 J^ P m 2 ^o + ^%)dx 2 . (3.14) 

By further assuming is compactly supported in either (—1,0) or (0, 1), we find that 
ip satisfies the equation (3.1) in a weak sense on (—1,0) and (0, 1). A standard bootstrap 
argument then shows that ip is in H k ((-1, 0)) (resp. H k ((0, 1))) for all k > when restricted 
to (—1,0) (resp. (0,1)), and hence is smooth when restricted to either interval. Since 
ip G Hq((—1, 1)), the jump conditions (3.2) and the boundary conditions (3.5) follow trivially. 
It remains to show that the jump conditions (3.3) and (3.4) hold. For this we take ip £ 
C^°((— 1, 1)) in (3.14). Integrating (3.14) by parts the terms of ipQ, we obtain 

f sMKiVih + m Wo + ki Wb - v%) + \B\ 2 m 2 4>'4>o - « ^ 2 (3.15) 
= mki v + r) + isi vmw + iei 2 b]^(o)^o(o) - a 2 jf 1 P (iei 2 ^o + V'Vo) 

Integrating further (3.15) by parts the terms of -?/> , using the fact that ip solves (3.1) on 
(— 1, 0) and (0, 1) and the jump conditions (3.2), we find that 

-mici 2 ^ + m + wu(o) + m<t - m V) + i^vi^o) 

= ([AW]+»MieiV(0))^b(0). (3.16) 

Since ip may be chosen arbitrarily, (3.16) implies the jump conditions (3.3)-(3.4) and we 
conclude our proof. □ 

Now we come to the heart part of this paper to clarify the sign of the infimum obtained 
in Proposition 3.1. It is crucial to represent the energy E[ip) in the following form 

Ety) = \t\ 2 E (il;) + 8E 1 y;), (3.17) 

where 

EoM = \ J\ W (V | 2 + ^) dx 2 - \g[p}\iP(0)\ 2 , (3.18) 

EiM = \ f M4|£| Vl* + I ^1 V + r\ 2 ) dx 2 . (3.19) 

Since the parameter s is positive but can be made to be small arbitrarily, the key point is to 
clarify the sign of the energy E (ip). Note that if without magnetic field this energy is non- 
positive, hence the pure fluid Rayleigh- Taylor problem is unstable. However, the presence of 
the magnetic field makes it have possibility to be positive. It is characterized by the critical 
magnetic value \B\ C and the critical frequency |£|f c (or |£|f c for the horizontal case), we recall 
their definitions (2.4)-(2.6). 
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Lemma 3.2. (i) The supremum (2.4) is achieved; 

(ii) For any fixed B so that \B\ < \B\ C , the infimum (2.5) is achieved. Moreover, |£|^ c is 
continuous and strictly decreasing as a function of \B\ within < \B\ < \B\ C and we have 

|f|f c -»0 as \B\->0, and |£|f c ->■ oo as \B\ -> \B\ C . (3.20) 

(Hi) For any fixed B so that \B\ < \B\ C , the supremum (2.6) is achieved. Moreover, |£|^ c 
is continuous and strictly decreasing as a function of \B\ within < \B\ < \B\ C and we have 

|£|f c ^oo as \B\^0, and ->■ as \B\ ->■ \B\ C . (3.21) 

Proof. To prove (i), we express the supremum (2.4) in an equivalent form 

i r 1 

= inf J \^'\ 2 dx 2 , (3.22) 



\B\l 

where ip e ifg((— 1, 1)) satisfying the constraint 

g[p}^ 2 (0) = 1. (3.23) 

It is obvious that the integral in (3.22) is bounded below. Let ip n be a minimizing sequence, 
then ip n is bounded in Hq and hence ip n — > ij) weakly in Hq, up to the extraction of a 
subsequence if necessary. Owing to the compact embedding Hq m> C° we have ip n (0) — > ip(0). 
By these convergences and the weak semi-continuity of the integral in (3.22), we obtain that 
if) is a minimizer of (3.22) and this proves (i). 

Now we prove (ii). To prove that the infimum (2.5) is achieved, we first rewrite (2.5) 
equivalently 

(lelfj 2 := mf PW, (3.24) 



where 



P(*P) = \B\ 2 J' W\ 2 dx 2 , (3.25) 

C = E H 2 ((-l,l)) g[p]^(0) - \B\ 2 £ W\ 2 dx 2 = 1 j . (3.26) 

Let ipn £ C be a minimizing sequence, we have from (3.24) that ip" is bounded in L 2 . This 
with Poincare's inequality implies that ip n is bounded in Hq. So we have tp n — > ip weakly in 
Hq, up to the extraction of a subsequence if necessary, and owing to the compact embeddings 
Hq ^ C 1 we have ip n — > ip strongly in H 1 and tp n (0) — > ip(0) as well. These convergences 
and the weak semi- continuity yield 

P(ip) < lim sup P($ n ) = inf P. (3.27) 

n— >oo C 

ip G C follows from the strong convergences and hence is a maximizer of (3.22). On 
the other hand, by the definitions (2.4) and (2.5), it is easy to have the conclusions of the 
behavior of |£|f c as B varies within < \B\ < \B\ C . 
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At last, we prove (Hi). We first rewrite the supremum (2.6) equivalently 

(|£|f c ) 2 :=supQW, (3.28) 



where 



Q(V) = g[ P mO) - \B\ 2 J* dx 2 , (3.29) 



(3.30) 



We first show that Q is bounded above on J 7 . Indeed, for any ip e J 7 , we have 



<C J^\d X2 (i) 2 )\dx 2 -\B\ 2 J \i)'\ 2 dx 



< C y ete 2 < C. (3.31) 

Let ^ n e J be a maximizing sequence, the second inequality in (3.31) implies that 

K j Wn\ 2 dx 2 <C J m 2 dx 2 -Q(ij n )<C. (3.32) 

Then ^ n is bounded in ifg, so we have ip n ip weakly in ifg, up to the extraction of a 
subsequence if necessary, and owing to the compact embeddings Hq y L 2 and Hq >• C° we 
have ip n ip strongly in L 2 and ^ n (0) — >■ ^(0) as well. These convergences and the weak 
semi-continuity of the integral in Q, we have that 

Q(tp) > lim sup Q (ip n ) = supQ- (3.33) 

n— >oo T 

ip G J 7 follows from the strong L 2 convergence and hence ip is a maximizer of (3.28). Finally, 
by the definitions (2.4) and (2.6), we conclude the behavior of |£|f c as B varies within 
< \B\ < \B\ C . The proof of Lemma 3.2 is completed. □ 

Now we can clarify the sign of the energy Eq(i/j). 

Lemma 3.3. We have the following four assertions: 

(i) If \B\ > \B\ C , then for any tp we have E (ip) > 0. Moreover, 

^oW>^yV| 2 -|5|^)l^1 2 + ^^^ 2 . (3.34) 

(ii) Fixed \B\ < \B\ C and |£| < |£|^ c; then for any vp we have E (ip) > 0. 

(Hi) Fixed \B\ < \B\ C and |£| > |£|^ c , then there exists ip such that Eo(ip) < 0. 
(iv) If there is ip such that E (ip) < 0, then \B\ < \B\ C , |£| < |£|f c and ip(0) ^ 0. 
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Proof. The first three assertions follows easily by the definitions (2.4) and (2.5). Indeed, 
for (i), we let B be so that \B\ > \B\ C . By the definition (2.4), we have 

\ B \l J W?dx 2 > g[p]^ 2 (0), for any ij) G H%. (3.35) 

Since the remaining term in Eq is nonnegative, we verify the assertion (i). To prove (ii) and 
(Hi), we fix any B so that \B\ < \B\ C . Again similarly by the definition (2.5), we have that 
if |e| < \Z\* C then 

f W(\eW? + m 2 ) dx 2 > g\Z\ 2 \pM0)\ 2 - (3.36) 

This proves the assertion (ii). Now if |£| > |£|f c , by the definition (2.6) there exists ip G Hq 
such that 

f W(\i\ 2 W\ 2 + WV) dx 2 < g\i\ 2 [p]\m\ 2 ■ (3.37) 

This prove (Hi). 

It remains to prove (iv). If E (ijj) < 0, then the assertions \B\ < \B\ C and |£| > |£]^ c 
follow from (i) and (ii). Since the integrals in (3.18) are all non- negative, then we must have 
ijj(0) 0. This proves (iv) and we conclude our lemma. □ 

By Lemma 3.3 we are able to show the sign of the infimum of E over A. We write 
a = a(s) to emphasize the dependence on s G (0, oo). 

Lemma 3.4. (i) // either \B\ > \B\ C or \B\ < \B\ C with |£| < |^|f c; then a(s) > for any 
s > 0. 

(ii) If both \B\ < \B\ C and \£\ > \£\y C are satisfied, then there exists s > depending on 
the quantities p±,[x±,g, \B\, |£| so that for s < s Q it holds that a(s) < 0. 

Proof. Since both E and J are homogeneous of degree 2 and J is positive definite, we may 
reduce to clarify the sign of the energy E(ip). 

To prove (i), observe that in this case we deduce from Lemma 3.3 that E (ip) > for any 
tp G A. Combining this with the fact E 1 (tp) > 0, we have E(ip) > for any s > 0. Hence 
taking the infimum we have fi(s) > 0. This proves (i). 

It remains to prove (ii). In this case we know from Lemma 3.3 that there exists ip such 
that E (t/j) < 0. Obviously, we have 

E$) = \Z\ 2 E (4>) + sE 1 $) < Eq$) + sC (3.38) 

for a constant C depending on p±,\i±,g, \B\, |£|. Then there exists So > depending on 
these parameters such that for s < sq it holds that E(ip) < 0. Hence the infimum a(s) < 
for s < sq. This proves (ii) and the proof of Lemma 3.4 is completed. □ 

Immediately, Lemma 3.4 proves the assertions (i) and (ii) of Theorem 2.1 by contradic- 
tion. To finish the proof of Theorem 2.1, it suffices to prove the assertion (Hi). So fixed 
\B\ < \B\ C and then |£| > |£|f c . We want to show that there is a fixed point such that A = s. 
To this end, we first study the behavior of a(s) as a function of s > 0. 
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Lemma 3.5. We have the following statements, 
(i) ac(s) is strictly increasing; 

W^e^ c 1 ((o,oo))nc°((o,oo)) ; 

(Hi) For any b > there exist two constants Cq,C\ > depending on the parameters 
P±, fi>±, 9, \B\, b so that 

a(s) < -C + sC x , for all |f | G [b, oo); (3.39) 

(iv) There exist constants C 2 > depending on p±,g and C3 > depending additionally 
on fj,±, \B\,£ so that 

a( S )>-C 2 \£\ + sC 3 . (3.40) 

Proof. Recall the energy decomposition (3.17) along with (3.18) and (3.19). It keeps the 
same form as in [5, Propostion 3.6], hence (i) and (ii) follow in the same way. 

To prove (Hi), fixed b > |£|^ c , by Lemma 3.3 (Hi) there exists ip b such that C = 
—E (ip b ) > 0. Then we have E(ip b ) < —C + sC\ for some C\ > and then (Hi) holds. 

Finally, we prove (iv). First observe that for any ip e A we have 

12 \ ltUf„l / / lelL/,|2 \ I I L/,/12 



-iei^Hi^(o)i 2 >-\t\g\p)[j m\ 2 dx 2 j u m 2 dx 2 j >-c a iei. (3.41) 

Since the other terms in the energy E is nonnegative, we have 

a(s) > -C 2 \C\ + s inf E x {$). (3.42) 

We denote by C3 this positive infimum, then (iv) follows and we conclude our lemma. □ 
By Lemma 3.4 and Lemma 3.5, we then define the open set 

S = ^((-00, 0)) C (0, 00). (3.43) 

Note that S is non-empty and allows us to define X(s) = ^J—a(s) for s G S. We state the 
existence of solutions to the modified problem (3.1)-(3.5) which we have already proved. 

Proposition 3.6. For each s G S there is a solution if> = ip s (^,x 2 ) with A = A(£, s) > to 
the problem (3.1) along with the jump and boundary conditions (3.2)-(3.5). Moreover, ip s , A 
are even in £. The solutions are smooth when restricted to (—1, 0) or (0, 1) with ip s (£, 0) 7^ 0. 

Proof. Let i/j s (£, x 2 ) be constructed in Proposition 3.1. Since s G S, we can write — a(£, s) = 
A 2 (£, s), then i[) s (^,x 2 ), A(£, s) solve the problem (3.1)— (3.5). The remaining assertions follow 
from Proposition 3.1 (ii) and Lemma 3.3 (iv). The proof of lemma is completed. □ 

Now we will make a fixed-point argument to find s G S such that s = A(|£|, s) to construct 
solutions to the original problem (1.48)-(1.52). 



Lemma 3.7. There exists a unique s G S so that A(|£|, s) = y/—fi(s) > and 

s = \(\£\,s). (3.44) 
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Proof. By Lemma 3.5, there exists s* > such that 

S = ^ 1 ((-oo,0)) = (0,s*). (3.45) 

We define A = y/—/J, on S and define the function $ : (0, s*) — >■ (0, oo) by 

$( S ) = S /A(|e|,s), (3.46) 

which is continuous and strictly increasing in s. Moreover, lim s ^o &( s ) = and lim^^ &(s) — 
+oo. Hence there is unique s G (0, s*) so that $(s) = 1, which gives (3.44). The proof of 
Lemma 3.7 is completed. □ 

In view of Proposition 3.6, Lemma 3.7 and Lemma 3.4, we conclude Theorem 2.1. □ 
3.2 Proof of Theorem 2.2 

The proof of Theorem 2.2 is similar to that of Theorem 2.1. The strategy is the same and 
the only difference is the energies defined when using the variational method to construct 
solutions. We define the two energies related to the problem (1.56)-(1.60) by 

E&) = \ f s^m 2 w\ 2 + 1 iei 2 ^ + r\ 2 ) + \Bmm 2 + \em 2 ) ^ 
—legmm 2 , (3.47) 

m = \ j\m 2 \M 2 + W\ 2 )dx 2) (3.48) 

which are both well-defined on the space H%((— 1, 1)). 
We define 

KW = \ £ VI 2 + W\ 2 ) dx 2 - \g[ P ]\m\ 2 - (3.49) 

As in the proof of Theorem 2.1, to prove Theorem 2.2 the first thing is to clarify the sign of 
E' (ip). We have the following lemma. 

Lemma 3.8. We have the following three assertions: 

(i) If \B\ > \B\ C , then for any ip we have E (tp) > 0. Moreover, 

E'M > \ J\\B\* ~ \B\l)W\ 2 + \B\^\^dx 2 . (3.50) 

(ii) Fixed \B\ < \B\ C and |£| > |£|f c , then for any ip we have Eo(ip) > 0. 

(Hi) Fixed \B\ < \B\ C and < |£| < |£|f c , then there exists ip such that E' (i[)) < 0. 

Proof. The proof follows by the definitions (2.4) and (2.6), which is similar to the that of 
Lemma 3.3. □ 

Once Lemma 3.8 is established, Theorem 2.2 follows similarly by the proof of Theorem 
2.1. □ 
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3.3 Proof of Theorem 2.3 

In this subsection we assume that \B\ < \B\ C to prove Theorem 2.3. Notice that in Lemma 3.7 
the fixed point s G S is unique, we may write uniquely A(|£|) within (|£|f c , oo) (resp. (0, |£|f c )) 
for B is vertical (resp. B is horizontal), while the corresponding solution constructed in 
Theorem 2.1 (iii) (resp. Theorem 2.2 (Hi)) is written by V^i- 

The continuity assertion follows in the same way as in [3 Proposition 3.9], it suffices to 
prove (2.7) and (2.8). To prove (2.7), we first derive the limit behavior. For this we take 
any |£| n G (|£|f c , oo) so that |£|„ — >■ |£|f c , then by Theorem 2.1 (iii) there exist functions 
V>|g| n G A so that 

-A 2 (|eW = i?(%|J<0. (3.51) 
Recalling the expressions of energies (3.17), (3.18) and (3.19), we have 

l£n| 2 #[p] ,2 , n \ \ B \ 2 f 1 nt- |2|/ |2 , I ;.// |2>, 



< A^(ieU) < ^Y^ k (0) " V 7 ^^11^1 + WJ ) dx ^ (3.52) 

We deduce from (3.51) and the fact ip\^\ n G .4. that is uniformly bounded in Hq as in 
the proof of Lemma 3.2. So, up to the extraction of a subsequence we have that 

i> weakly in ifo((-l, x )) and strongly in H%((-1, 1)) (3.53) 

as well as ^^(O) — > ip(0). So taking superior limit in (3.52) as n — > oo along the subsequence, 
we have 

< hmsu P A 2 (KU) < (|flS!gM ^(o) - H! [\(\Z0 2 m 2 + Wf)dx2 < 0. (3.54) 

The last inequality above comes from Lemma 3.3 (ii). Since (3.54) holds for any such 
extracted subsequence, we deduce that lim^oo A 2 (|£| n ) = for the original sequence |£|„ as 
well. This proves lim^^is A(|£|) = 0. 

Now we turn to prove the boundedness of A. Using the fact that — A(|£|) 2 = E(ip\^\) and 
the expressions of energies (3.17), (3.18) and (3.19), we have 

A(iei) 2 < ^l£| 2 <7[p]|%|(0)| 2 - ~ £ \B\ 2 m 2 W^\ 2 + |^|| 2 ) dx 2 . (3.55) 
Notice additionally that G A, we have 

= \j\m 2 \M 2 + W)d*2 = 1. (3.56) 
Hence from (3.55)-(3.56) we can bound A by 

i i 
A(KI) 2 <y0^=(\B\*J o \em 2 dx,y Qf P+ \e\^\\ 2 dxi 

|2/|t|2|„/./ |2 i I „ /.// |2\ 



-1 



\Bmw m r + \n\ndx2 



^Jl^^-Mi- (3 - 57) 
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Hence, (2.7) follows. 

Finally, to prove (2.8), the limit lim^i^^s A(|£|) = can be proven in the same way as 
lim^i^iB A(|£|) = 0. It remains to prove lim^^o A(|£|) = and this follows from (3.55)- 
(3.56) again by 

The proof of Theorem 2.3 is completed. □ 
3.4 Proof of Theorem 2.4 

In this subsection we will do the energy estimates to prove Theorem 2.4. In this subsection 
and also in Theorem 2.4, L 2 , H l denote the usual LP and Sobolev spaces on Q, and their 
norms are denoted by || -\\l 2 , || • Wh 1 respectively. We denote C be generic constants depending 
only on the physical coefficients and B. 

First, we prove (i) and let B = (0,1?) for the moment. To prove (2.10), differentiating 
the second equation in (1.35) with respect to time t and eliminating the rj term by using the 
first equation, we obtain 

pd tt v + Vd t q - fiAdtV - \B\ 2 d 2 22 v = 0, _ _ 

divy = div<9ft> = 0, 

along with the jump and boundary conditions 

M = IdM = 0, {-n(DdtV + Dd t v T ) + d t qlje 2 = g[p}v 2 e 2 + \B\ 2 {d 2 vl (3.60) 

v_(t,X!, -1) = v+(t,xi, 1) = d t v_(t,x 1 , -1) = d t v + (t,x 1 , 1) = 0. (3.61) 
We regard the system (3.59)-(3.61) as one for dtv, dtq, v and we have the energy identity: 

Lemma 3.9. Let d t v,d t q,v solve (3.59)-(3.61), then we have 

^^pI^^ + isH^i'^-^mMx^o)! 2 ^ 

+ / fi\Dd t v + Dd t v T \ 2 dx = 0. (3.62) 
Jn 

Proof. Multiplying the equations (3.59)i by d t v and integrate over f2, after integrating by 
parts respectively in Q + and f2_ and using the jump and boundary conditions (3.60)-(3.61), 
by (3.59) 2 , we obtain (3.62). □ 

Lemma 3.10. Suppose \B\ > \B\ C , then there exists constant C^ 1 > such that 

\B\ 2 \d 2 v\ 2 dx - [ ^]Mxi,0)| 2 da;i>C- 1 (||7;||i2 + ||^||i2). (3.63) 



J 

Jn 
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Proof. We represent the integral above as 

/ \B\ 2 \d 2 v\ 2 dx - [ sMhO^o)! 2 ^! 

Jn Jr 

= j£(/ \B\ 2 \d 2 v(x 1 ,x 2 )\ 2 dx 2 -g[p]\v 2 (x 1 ,0)\^j d Xl . (3.64) 

By the definition (2.4) of the critical magnetic number \B\ C , we have that for any x\ G R, 



l 

2|Q „./„ ™ M2 j„ „r„lL. n\|2 



1 5 1 | d 2 v (xi , x 2 ) r ^2 - 9 H I V 2 (Xi , 0) | ' 

1 

> y (\B\ 2 \d 2 v 1 (x 1 ,x 2 )\ 2 + (\B\ 2 -\B\ 2 c )\d 2 v 2 (x 1 ,x 2 )\ 2 ) dx 2 . (3.65) 

Substituting (3.65) into (3.64), we have 

/ \B\ 2 \d 2 v\ 2 dx - [ ^[p]h(xi,0)| 2 tfai>C- 1 ||^u||i a . (3.66) 
Jn Jr 

Hence, (3.63) follows from (3.66) and Poincare inequality in the slab, since the boundary 
condition (3.61). The proof of Lemma 3.10 is completed. □ 

Now, we define 

£(dtv,v;t):= [ p\d t v\ 2 + \B\ 2 \d 2 v\ 2 dx - [ g\p]\v 2 (x u 0)\ 2 dx u (3.67) 
Jn Jr 

V(d t v; t):= I n \Dd t v + Dd t v T \ 2 dx, (3.68) 
Jn 

then by Lemma 3.10, we have 

C-WldtvWl? + \\v\\ 2 L2 + \\d 2 v\\ 2 L2 ) < S(dtv,v; t) < C(\\d t v\\ 2 L2 + \\v\\ 2 L2 + \\d 2 v\\ 2 L2 ), (3.69) 
and by Korn's inequality in the slab, since the boundary condition (3.61), we have 

C-lSH&i < V{d t v-t) < C\\d t vf m , (3.70) 
Hence, integrating (3.62) directly in time, by (3.69)-(3.70), we obtain 



\\dtv{t)\\l* + \\v{t)f L2 + UMOIIi* + I \\dtv(s)\\ 2 m ds 
< C{\\d t vm\h + IK0)||i 2 + 11^(0)11^)- (3-71) 

This proves (2.10). 

To prove (2.11), we notice that we have the boundary condition for rj, 

r j -(t,x 1 ,-l) = r j+ (t,x 1 ,l)=0, (3.72) 
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which is deduced from the first equation in (1.35) and the initial assumptions on rf. Then 
v,q,t] satisfy the system (3.59)-(3.61) by replacing d t v,d t q,v correspondingly. Hence, the 
arguments for proving (2.10) also lead to (2.11). 

Now we turn to prove (ii) and hence let B = (-5,0). The proof is similar to that of (i). 
Similarly, to prove (2.12) we have the following system 



pd tt v + V8 t q - pAd t v - \B\ 2 d 2 n v = 0, 
divv = div<9 t i> = 0. 



(3.72) 



along with the jump and boundary conditions 

H = {dtvj = 0, \-n{DdtV + Dd t v T ) + d t qlje 2 = g\p]v 2 e 2 , (3.73) 

v-(t,xx, -1) = v+(t,xi, 1) = d t v-(t,xi, -1) = d t v+(t,x 1 , 1) = 0. (3.74) 
We have the energy identity for the system (3.72)-(3.74): 

Lemma 3.11. Letv,q solve (3.72)-(3.74), then we have 



^(/ Pl ^ |2 + |B|2| ^ |2 ^~I^ P]h(;ri,0)|2rfXl 

+ / /i\Dd t v + Dd t v T \ 2 dx = 0. (3.75) 
Jn 

Proof. Multiplying the equations (3.72) x by v and integrate over Q, after integrating by 
parts respectively in Q + and f2_ and using the jump and boundary conditions (3.73)-(3.74), 
by (3. 72) 2 , we obtain (3.75). □ 

Lemma 3.12. Suppose \B\ > \B\ C , then there exists constant C^ 1 > such that 

[ \B\ 2 \d lV \ 2 dx - [ g[p]\v 2 (x 1 ,0)\ 2 dx 1 >C- 1 (\\d 1 v 1 \\ 2 L2 + \\v 2 \\ 2 Hl ). (3.76) 
Jn Jr 

Proof. By d±vi = —d 2 v 2) we can represent the integral in (3.76) as 

f \B\ 2 \d lV \ 2 dx - [ gipWv^x.M 2 dx, 
Jn Jr 

= [ \B\ 2 (\d 2 v 2 \ 2 + \d lV2 \ 2 )dx- [ g[p]\v 2 (x 1 ,0)\ 2 dx 1 
Jn Jr 

= jT \B\ 2 (\d 2 v 2 \ 2 + \d 1 v 2 \ 2 )dx 2 -g[p}\v 2 (x 1 ,0)\ 2 ^ dx ± . (3.77) 
By the definition (2.4) of \B\ C , we have 

J \B\ 2 \d 2 v 2 ( Xl ,x 2 )\ 2 dx 2 - g[p}\v 2 ( Xl ,0)\ 2 

> j (\B\ 2 -\B\ 2 c )\d 2 v 2 (x u x 2 )\ 2 dx 2 . (3.78) 
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Substituting (3.78) into (3.77), we have 

f \B\ 2 \d lV \ 2 dx - [ g^ix^o^dx^c-WidmWh + M 2 ^). (3.79) 

Jn Jr 

Hence, (3.76) follows from (3.79) and Poincare inequality in the slab, since the boundary 
condition (3.74). The proof of Lemma 3.12 is completed. □ 

Now, we define 



£'(d t v,v;t) := [ p\d t v\ 2 + \B\ 2 \d lV \ 2 dx - [ g[p]\v 2 (x u 0)| 2 dx u (3 
Jn Jr 



80) 



then by Lemma 3.12, we have 

C'WldtvWh + \\d lVl \\l 2 + \\v 2 \\ 2 m ) < £'(d t v,v;t) < C(\\d t v\\ 2 L2 + \\d lVl \\ 2 L , + \\v 2 \\ 2 m ).(3.81) 
Hence, integrating (3.75) directly in time, by (3.81) and Korn's inequality, we obtain 



\\dMt)\\h + \\dMt)\\h + IMOIIk + jT \\d t v(s)\\ 2 Hl ds 

< C(\\dM0)\\h + ll^i(0)||| 2 + \\v 2 (0)\\ 2 Hl ). (3.82) 

This proves (2.12). 

To prove (2.13), by the initial condition and the first and third equations in (1.35), 

dlvr] = 0. (3.83) 

Hence, the situation is the same to the proof of (2.12) and hence along the similar lines we 
obtain (2.13). The proof of Theorem 2.4 is completed. □ 
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